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I.  INTRODUCTION 

This  report  describes  the  work  done  on  Contract  No.  N00123-76-C- 
1371  for  the  period  1  November  1976  to  31  January  1977. 

The  overall  program  is  divided  into  three  areas.  These  are 
1)  basic  scattering  code  development,  2)  reflector  antenna  code  develop¬ 
ment  and  3)  basic  theoretical  studies  to  support  the  first  two  areas. 

The  following  sections  describe  the  progress  made  during  the  second 
quarter  in  each  of  the  three  areas  mentioned  above. 


II.  PROGRAM  SCOPE 

The  scope  of  the  work  under  Contract  No.  N001 23-76-C-l 371  is  to 
develop  the  necessary  theory,  algorithms  and  computer  codes  for  simu¬ 
lating  antennas  at  UHF  and  above  in  a  complex  ship  environment.  A 
milestone  chart  for  the  total  program,  which  extends  over  a  three  year 
period,  is  shown  in  Table  I.  A  more  detailed  breakdown  of  the  effort 
planned  for  the  first  year  is  shown  in  Table  II.  The  following  sections 
describe  the  progress  made  during  the  second  quarter  of  the  program  in 
Table  II. 


III.  BASIC  SCATTERING  CODE  DEVELOPMENT 

The  purpose  of  this  section  is  to  describe  the  present  status  of 
the  basic  scattering  code  development  for  the  analysis  of  antennas  in 
a  complex  shipboard  environment.  The  Geometrical  Theory  of  Diffraction 
(GTD)  is  being  used  to  develop  algorithms  to  solve  for  the  scattering 
from  basic  plate  and  cylinder  structures.  These  simple  components  can 
then  be  combined  to  form  box-like  structures  and  finite  elliptic  cylin¬ 
ders  that  can  represent  various  component  structures  of  a  ship.  The 
algorithms  are  being  implemented  into  a  user-oriented  computer  code. 

In  this  period  the  multiple  plate  scattering  code  has  been  improved 
and  tested.  This  multiplate  code  allows  a  number  of  N  sided  flat  plates 
to  be  placed  in  proximity  to  each  other  such  that  they  can  form  arbitrary 
box-shaped  structures.  Some  improvements  have  been  made  in  the  algorithms 
to  help  the  accuracy  and  efficiency  of  the  computer  codes.  For  example, 
the  algorithm  for  finding  reflected  fields  from  flat  surfaces  has  been 
improved  by  a  more  efficient  image  theory  solution.  Also,  the  algorithm 
used  to  find  the  incident  field  needed  for  the  slope  diffracted  field 
has  been  improved  such  that  it  will  allow  use  of  a  more  general  source. 

The  validity  of  the  multiple  plate  scattering  code  is  being  tested 
by  comparing  the  results  against  experimental  measurements.  The  first 
comparisons  were  made  with  a  dipole  antenna  near  two  plates  forming  a 
convex  edge  (roof-like  structure)  as  illustrated  in  the  insert  of  Fig.  1. 
The  GTD  and  measured  results  are  compared  in  Figs.  1-5  for  various  pat¬ 
tern  cuts  and  polarizations  as  illustrated  in  the  figures.  The  agreement 
is  very  good  in  all  cases. 
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RELATIVE  POWER  (dB) 


One  of  the  reasons  for  this  excellent  agreement  in  all  tne  pattern 
cuts  is  the  use  of  newly  developed  empirical  corner  and  slope-corner 
diffraction  coefficients  (vertex  diffraction).  To  show  how  important 
these  terms  are  in  finding  volumetric  patterns.  Figs.  6  and  7  show  GTD 
results  for  the  above  roof-shaped  structure  without  the  corner  diffrac¬ 
tion  terms  compared  against  measured  results.  The  discontinuities  shown 
in  the  patterns  result  from  the  regular  edge  diffraction  and  slope  dif¬ 
fraction  terms  "falling  off"  at  the  f i n i .  _  limits  of  the  edges  of  the 
plates.  The  empirical  corner  and  slope-corner  diffraction  terms  smooth 
out  these  discontinuities  as  was  previously  shown  for  these  two  cases 
in  Figs.  2  and  3,  respectively. 

It  must  be  cautioned  that  these  newly  developed  empirical  corner 
diffraction  terms  have  not  been  put  on  a  firm  rigorous  foundation  as 
yet.  This  should  be  done  in  order  to  have  complete  confidence  in  the 
results  for  all  possible  applications.  The  primary  justification  up  to 
this  point  is  the  experimental  agreement.  It  is  felt  that  this  is  suffi¬ 
cient  evidence  to  continue  their  use  at  the  present  time. 


MEASURED 


Fig.  6.  Comparison  of  the  measured  and  calculated  (without 
corner  diffraction)  E.  radiation  pattern  (refer 
to  Fig.  2). 
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Fig.  7. 
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Comparison  of  measured  and  calculated  (without 
corner  diffraction)  E,  radiation  pattern 
(refer  to  Fig.  3). 
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The  scattering  code  also  has  been  tested  on  a  box-shape  as  illustrated 
in  the  insert  of  Fig.  8.  The  scattering  code  result  without  the  corner 
diffraction  terms  is  compared  with  measured  results  in  Fig.  8a.  Note 
that  the  empirical  corner  diffraction  terms  are  again  needed  to  obtain 
a  continuous  field  pattern.  This  is  shown  in  Fig.  8b  which  compares  the 
scattering  code  result  with  the  empirical  corner  diffraction  coefficient 
included  against  measured  results.  Comparison  between  the  GTD  with 
corner  diffraction  and  measured  results  are  shown  for  two  other  cases  in 
Figs.  9  and  10.  The  good  agreement  obtained  again  confirms  the  validity 
of  the  theoretical  results. 

The  basic  scattering  code  for  a  finite  elliptic  cylinder  also  has 
been  investigated  during  this  period.  In  particular,  the  diffraction 
mechanisms  used  to  model  the  junction  between  the  flat  plate  endcaps  and 
the  curved  elliptic  cylinder  surface  have  been  studied.  A  curved  edge 
dyadic  diffraction  coefficient  developed  by  Kouyoumjian  and  PathakpJ  is 
used  to  find  the  fields  scattered  from  the  edges.  In  general,  there  are 
four  points  of  diffraction  around  one  of  these  elliptic  rims  for  a  given 
observation  direction.  Previously,  a  search  technique  was  used  to  locate 
these  points.  However,  this  proved  to  be  a  very  time  consuming  means  of 
locating  the  points  of  diffraction.  A  faster  method  has  been  developed 
for  finding  these  points  by  numerically  solving  for  the  roots  of  an 
eighth  order  polynomial  with  coefficients  that  are  a  function  of  the 
elliptic  rim  parameters  and  observation  direction.  This  is  essentially 
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Comparison  of  measured  and  calculated  (without  corner 
diffraction)  Ee  radia’on  pattern  for  a  dipole  near 
a  box  in  the  indicated  plane. 
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an  exact  solution  of  the  diffraction  points  within  the  accuracy  of  the 
algorithms  for  the  numerical  polynomial  solution  routines.  The  speed 
of  finding  the  diffraction  points  has  been  increased  by  a  factor  of  3.75 
and  consequently  the  overall  program  speed  has  been  improved  by  a  factor 
of  approximately  2.  The  computer  storage  needed  for  the  program  also 
has  been  reduced  in  the  process.  Further  improvements  in  the  scattering 
code  are  under  investigation. 

In  the  next  period,  the  scattering  code  for  the  plates  and  box¬ 
shaped  structures  will  be  completed.  It  will  be  delivered  to  NELC 
early  in  the  fourth  quarter.  This  early  delivery  of  the  first  version 
of  the  scattering  code  will  allow  for  inputs  from  NELC  on  improvements 
of  the  usefulness  of  subsequent  scattering  codes- 


The  scattering  code  for  the  isolated  finite  elliptic  cylinder  will 
be  improved  and  tested  as  has  been  done  for  the  plate-box  code.  When 
this  is  completed  the  scattering  code  for  the  cylinder  will  be  combined 
with  the  plate-box  code,  but,  it  will  not  contain  box-cylinder  inter¬ 
actions  in  this  first  combined  code.  This  combined  scattering  code  will 
be  delivered  sometime  at  the  end  of  the  fourth  quarter  along  with  a  users 
manual . 
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Fig.  10.  Comparison  of  measured  and  calculation  radiation 
pattern,  a)  E0  and  :■)  E(!). 


This  final  code  for  the  first  year  will  allow  field  pattern  com¬ 
putations  for  isolated  plates,  boxes,  finite  elliptic  cylinders,  and  a 
psuedo  combination  of  both  without  box-cylinder  interactions.  The  theory 
for  the  coupled  solution  that  will  include  the  dominate  interactions  for 
the  box-cylinder  far  field  patterns  will  be  under  investigation  for  the 
remainder  of  this  contract  year. 


IV.  REFLECTOR  ANTENNA  CODE  DEVELOPMENT 

The  purpose  of  the  present  effort  is  to  develop  a  user-oriented 
computer  program  package  by  which  the  far  field  pattern  of  a  typical 
Navy  reflector  antenna  can  be  calculated.  Feed  blockage  and  scattering 
effects  are  not  included  in  this  phase.  These  will  be  included  later 
(see  Table  I). 

The  theoretical  approach  -or  computing  the  far  field  pattern  of  the 
general  reflector  is  based  on  a  combination  of  the  Geometrical  Theory  of 
Diffraction  (GTD)  and  Aperture  integration  (AI)  techniques.  AI  will  be 
used  to  compute  the  main  beam  and  near  sidelobes;  GTD  will  be  used  to 
compute  the  wide-angle  sidelobes  and  the  backlobes.  To  implement  the 
computer  algorithms  based  on  these  theories,  efficient  ways  are  needed 
to  handle  calculations  involving  the  feed  pattern  and  the  reflector 
geometry.  It  is  planned  to  treat  the  geometry  of  the  reflector  rim  by 
using  best-fit  elliptical  or  linear  segments. 

Since  the  majority  of  Navy  reflector  antennas  have  parabolic  sur¬ 
faces,  only  the  class  of  parabolic  surfaces  will  be  implemented  in  the 
computer  code.  The  code  for  the  reflector  geometry  will  be  flexible 
enough  to  include  off-set  fed  reflectors  and  general  reflector  rim  shapes 
such  as  elliptical  and  rectangular  with  chopped  corners. 

An  efficient  way  to  handle  AI  is  necessary  because  of  its  time- 
consuming  nature  for  electrically  large  reflectors.  Fortunately,  AI  is 
needed  only  for  the  main  beam  and  near  sidelobe  regions  of  the  pattern. 
GTD  is  very  efficient  for  calculating  most  of  the  pattern.  Two  techni¬ 
ques  are  being  implemented  to  provide  efficient  AI:  the  polynomial-fit 
method  for  the  feed  pattern  and  the  use  of  overlapping  rectangular  sub¬ 
apertures  for  the  two-dimensional  numerical  integration  required  to  cal¬ 
culate  the  far  field  pattern. 

The  polynomial -fi t  method  is  used  to  obtain  an  analytic  fit  for 
each  measured  pattern  cut  of  the  volumetric  feed  pattern.  This  method 
provides  a  computationally  efficient  way  of  calculating  the  aperture 
field  without  requiring  large  amounts  of  computer  storage  for  the  mea¬ 
sured  feed  pattern.  Only  relatively  few  coefficients  need  to  be  stored 
for  essentially  complete  feed  pattern  information.  Furthermore,  the 
polynomial-fi t  method  has  the  advantages  of  flexibility  and  simplicity 
for  general  feed  patterns.  No  cut-and-try  procedures  are  needed;  the 
polynomial  coefficients  can  be  computed  automatically  from  the  measured 
feed  pattern  input. 
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As  reported  in  the  previous  Quarterly  report,  computer  subroutines 
were  developed  to  provide  a  polynomial  curve-fit  for  measured  feed  pat¬ 
terns.  This  method  was  demonstrated  for  typical  sum  and  difference  pat¬ 
terns.  Further  improvement  was  made  tin's  quarter  for  fitting  measured 
feed  patterns  in  the  spillover  region,  i.e.,  where  the  feed  radiation 
is  not  reflected  by  the  reflector.  At  the  low  field  levels  encountered 
in  the  spillover  region,  a  piecewise  linear  approximation  was  found  to 
be  a  better  representation  for  this  part  of  the  feed  pattern.  Titus  the 
computer  subroutines  were  modified  to  use  a  combination  of  polynomial  - 
fit  for  the  major  part  of  the  feed  pattern  and  a  piecewise  linear  fit 
for  most  of  the  spillover  region  of  the  pattern. 

The  polynomial -fi t  method  is  based  on  approximating  a  measured 
feed  pattern  by 


(1)  f(b)  *  I  B 

t=l 


yl 


Qo  ”<l > 
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The  polynomial  fit  is  obtained  for  4^  <p  «•  ^  as  shown  in  Fig.  11. 

The  feed  pattern  may  be  a  symmetrical  sum  pattern  (even  terms  only),  a 
symmetrical  difference  pattern  (odd  terms  only),  or  a  general  unsymmetric 
pattern  (both  even  and  odd  terms).  Thus,  in  Eq.  (1) 
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even  terms  only 

odd  terms  only 

both  odd  and  even  terms 


For  a  symmetrical  sum  pattern  (even  terms  only)  and  with  the  usual 
value  of  ii/ 1  =  0: 


(3) 


f(,)  =  B,  *  62(L) 


\2N-2 


where  N  is  the  number  of  measured  feed  pattern  values  that  are  fitted. 
Thus  substituting  the  N  known  feed  values  f(bxi )  into  Eq.  (3)  gives  N 
equations  and  N  unknown  coefficients: 


(4) 


f (0)  =  Bj 


f(*x2>  =  b1+b2  (j?)  *  •••  *  bn(^) 
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By  using  a  simultaneous  equation  subroutine  on  Eqs.  (4),  the  coefficients 
Be  are  solved.  One  of  the  computer  subroutines  performs  the  operation  of 
reading  in  measured  feed  pattern  values  then  solving  for  the  coefficients 
B;  and  storing  them.  Another  subroutine  can  then  be  called  whenever 
needed  to  calculate  the  feed  pattern  from  Eq.  (3). 

For  regions  of  the  feed  pattern  where  there  is  ripple  or  low 
values  a  piece-wise  linear  approximation  is  used.  One  subroutine  calcu¬ 
lates  the  slope  of  each  section  as  shown  in  Fig.  11  from 

(5)  Afl  _  f^i+l)  ~  f(»j) 


Another  subroutine  can  then  be  called  when  needed  to  calculate  the  feed 
pattern  in  the  piece-wise  linear  region  from 

(6)  fU)  =  fU  ,•)  +  £  (*-*<) 

1  1 

for  t|i.j  _  4.  i'-j  +  i  • 

A  result  using  the  combination  of  polynomial -fit  and  piece-wise 
linear  approximation  on  a  symmetrical  sum  pattern  (even  terms  only)  is 
shown  in  Fig.  12.  The  polynomial  fit  was  used  for  0  ^  ^  <_  85°,  and  the 
piece-wise  linear  approximation  was  used  for  85°  <  c  <_  180°. 

Although  GTD  can  be  used  to  calculate  most  of  the  far  field  pattern, 
the  main  beam  and  near  sidelobes  of  the  reflector  antenna  will  be  calcu¬ 
lated  by  the  aperture  integration  approach.  Since  the  required  integra¬ 
tion  is  usually  very  cumbersome  for  electrically  large  antennas,  an 
approach  based  on  using  overlapping  rectangular  subapertures  is  being 
used.  Each  subaperture  can  be  electrically  large,  i.e.,  several  wave¬ 
lengths  in  size.  This  will  minimize  the  number  of  subapertures  required. 

Preliminary  results  have  been  obtained  using  this  aperture  integra¬ 
tion  approach  for  a  line  source  with  an  exactly  known  far  field  pattern 
as  a  check  case.  It  was  found  that  the  use  of  overlapping  rectangular 
subapertures  each  with  a  triangular  subaperture  distribution  gives 
excellent  accuracy  with  relatively  few  aperture  field  samples.  The 
aperture  distribution  for  the  line  source  used  in  the  check  case  is 
shown  in  Fig.  13.  A  spacing  of  2\  was  used  between  "sample"  points  along 
the  20>,  aperture.  As  seen  from  Fig.  13  the  overlapping  subapertures 
(triangular  distribution)  give  a  good  piece-wise  linear  approximation  to 
the  20>.  aperture  distribution.  The  use  of  non-overlapping  subapertures 
each  with  a  uniform  distribution  also  is  shown  in  Fig.  13;  this  technique 
gives  a  "step"  approximation  and  thus  requires  considerably  more  "sample" 
points  than  the  use  of  the  overlapping  subaperture  to  achieve  the  same 
accuracy.  The  exact  far  field  pattern  is  shown  in  Fig.  14  along  with 


SUBAPERTURES  WITH 
UNIFORM  DISTRIBUTION 


Fig.  13.  Exact  aperture  distribution  and  approximation 
to  it  by  subaperture  techniques. 
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patterns  calculated  from  the  ovtjr  loop  ■  *kj  subapert  u>es  (triangular  dis¬ 
tribution)  and  non-overlapping  subapertures  (uniform  distribution] . 

The  result  using  "array"  sampling  is  also  shown  in  Fig.  14.  The  rec¬ 
tangular  subapertures  with  triangular  distribution  appear  to  be  superior 
to  other  aperture  integration  methods  such  as  "array"  sampling  of  the 
aperture  field  or  rectanqular  snbapertures  with  uniform  distribution' 
and  will  be  incorporated  into  our  general  reflector  code. 


V.  THEORETICAL  STUDIES 

Edge  diffraction  is  part  of  the  foundation  of  the  Geometrical 
Theory  of  Diffraction  (GTD).  A  uniform  solution  for  the  diffraction  of 
an  electromagnetic  wave  by  a  perfectly  conducting  wedge  has  been  obtained 
within  the  format  of  the  GTD.  If  the  field  incident  on  the  edge  does 
not  have  a  rapid  spatial  variation,  the  diffracted  field  is  directly 
proportional  to  the  incident  field  at  the  edge  and  it  can  be  calculated 
using  the  dyadic  diffraction  coefficient  given  by  Kouyoumjian  and  Pathak 
[1]. 

In  many  practical  problems,  however,  the  incident  field  has  rapid 
spatial  variation  at  a  diffracting  edge.  If  the  edge  diffracted  field 
is  calculated  in  the  usual  way[l]  there  will  be  a  discontinuity  in  the 
spatial  derivatives  of  the  total  field,  i.e.,  the  calculated  pattern 
will  be  continuous  but  have  a  "kink"  in  it.  Some  examples  were  given 
in  the  previous  report[2]. 

Thus,  the  theoretical  studies  on  this  contract  to  date  have  con¬ 
centrated  on  the  development  of  an  accurate  yet  simple  slope-diffraction 
coefficient.  Some  details  were  given  in  the  previous  report[2].  This 
study  is  essentially  completed  and  full  details  will  be  presented  in  a 
technical  report  that  is  in  preparation.  Numerical  examples  will  be 
included  to  illustrate  the  usefulness  of  slope  diffraction. 
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